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Abstract 

Given a control region Q on a compact Riemannian manifold M, we consider the 
heat equation with a source term g localized in Q. It is known that any initial data 
in L 2 (M) can be steered to in an arbitrarily small time T by applying a suitable 
control g in L 2 ([0, T] x fi), and, as T tends to 0, the norm of g grows like exp(C/T) 
times the norm of the data. We investigate how C depends on the geometry of Q. 
We prove C > d 2 /4 where d is the largest distance of a point in M from £2. When M 
is a segment of length L controlled at one end, we prove C < a*L 2 for some a* < 2. 
Moreover, this bound implies C < a^L 2 ^ where Lq is the length of the longest 
generalized geodesic in M which does not intersect The control transmutation 
method used in proving this last result is of a broader interest. 



1 The problem 



Let (M, g) be a smooth connected compact n-dimensional Riemannian man- 
ifold with metric g and boundary dM. When dM ^ 0, M denotes the in- 
terior and M = M U dM. Let dist : M — > R + denote the distance func- 
tion. Let A denote the (negative) Dirichlet Laplacian on L 2 (M) with domain 
D(A) = Hl{M)r\H 2 (M). 
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Consider a positive control time T, and an open control region f2. Let 1]o,t[xQ 
denote the characteristic function of the space-time control region ]0,T[xf2. 
The heat equation on M is said to be null- controllable (or exactly controllable 
to zero) in time T by interior controls on Q if for all uq G L 2 (M) there is a 
control function g E L 2 (Rx M) such that the solution u G C°([0, oo), L 2 {M)) 
of the mixed Dirichlet-Cauchy problem: 

d t u-Au= 1 g in ]0, TlxM, u = on ]0, TlxdM, (1) 

]0,T[xSl 

with Cauchy data u — u at t — 0, satisfies w = at t = T. For a survey on 
this problem prior to 1978 we refer to [Rus78]. For a recent update, we refer to 
[ZuaOl]. Lebeau and Robbiano have proved (in [LR95] using local Carleman 
estimates) that there is a continuous linear operator S : L 2 (M) — > Cq°(M. x M) 
such that g = Su yields the null-controllability of the heat equation on M in 
time T by interior controls on f2. 

The most striking feature of this result is that we may control the heat in 
arbitrarily small time whatever geometry the control region has. In this paper 
we address the following question: How does the geometry of the control region 
influence the cost of controlling the heat to zero in small time ? 

Now, we shall formulate this question more precisely and give references. 

Definition 1.1. For all control time T and all control region Q, the null- 
controllability cost for the heat equation on M is the best constant, denoted 
Ct,ci, in the estimate: 

IMU 2 (RxM) < Ct,Q.\\uo\\l^(M) 

for all initial data uo and control g solving the null-controllability problem 
described above. 

By duality (cf. [DR77]), C T ,n is also the best constant in the observation 
inequality for the homogeneous heat semigroup t i— > e tA : 

Vw e L 2 (M), ||e TA M ||L2(M) < C T ,n\\e tA uo\\ L 2(( 0iT)xQ ) . 

Lebeau and Robbiano's result implies the finiteness of the null-controllability 
cost for the heat equation on M for any control time and any control region. 
Emanuilov extended this result to more general parabolic operators in [Ema95] 
using global Carleman estimates with singular weights. When (M , g) is an open 
set in Euclidean space, this method was used by Fernandez-Cara and Zuazua 
in [FCZ00] to obtain the optimal time dependence of the null-controllability 
cost for small time, i.e.: 

< sup p 2 /A < lim inf T In C T a < hm sup T\nC T n< +oo (2) 

B p cM\n t^o 
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where the supremum is taken over balls B p of radius p. The lower bound 
is stated in section 4.1 of [ZuaOl] and it is based on the construction of a 
"very singular solution of the heat equation in (0, +00) x M. n " used in the 
proof of Theorem 6.2 in [FCZ00]. Note that the method used in theorem 1 
of [LR95] seems to fall short of the optimal time dependence. Actually, using 
the improved version of proposition 1 in [LR95] presented as proposition 2 in 
[LZ98], we have only been able to prove that hmsup T _ +0 T 7 \nC T ,n is finite for 
all 7 > 1. 

Indeed Seidman had already asked how violent fast controls are, and his first 
answer concerned heat null-controllability from a boundary region T C dM. 
In [Sei84], under the condition that the wave equation on M is exactly con- 
trollable by controls in T in time L, he computes an explicit positive value 
(3 such that lim sup T ^ T In Ct,t < PL 2 (we give more explanations on this 
geometric upper bound in section 2 after theorem 2.3). The positivity of 
lim inf T ^ T In Ct,t when M is an interval was subsequently proved by Giiichal 
in [Gui85], ensuring the optimality of Seidman's result with respect to the time 
dependence. Later, Seidman also addressed finite dimensional linear systems 
as well as the Schrodinger and plate equations (cf. the companion paper [Mil03] 
for more details and references). 



2 The results 

2.1 Lower bound 

Our first result, proved in section 3, generalizes and improves on the geometric 
lower bound of Fernandez- Cara and Zuazua: 

Theorem 2.1. The null- controllability cost of the heat equation for small time 
(cf. definition 1.1) satisfies the following geometric lower bound: 

lim inf Tin C Tn > sup dist(y, H) 2/ 4 (3) 

As put in [ZuaOl], such a lower bound follows from the construction of a 
"very singular solution of the heat equation". Our construction underscores 
that only a large but finite number of modes is needed. For a short control 
time T > 0, we consider a Dirac mass as far from Q as possible, we smooth 
it out by applying the homogeneous heat semigroup for a very short time 
(eT with small e) and truncating very large frequencies (larger than (eT)~ r ), 
and finally we take it as initial data in (1). The proof relies on Varadhan's 
formula for the heat kernel in small time (cf. [Var67]), which requires very low 
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smoothness assumptions as proved in [Nor97]. 

We believe that there is no solution of the heat equation which is more singular 
than the heat kernel and therefore conjecture that this lower bound is also an 
upper bound, i.e. hinTlnC-m = supdist(y, f2) 2 /4. 

2.2 The segment controlled at one end 

Our second result, proved in section 4, concerns the most simple heat null- 
controllability problem: the heat equation on a segment controlled at one end 
through a Dirichlet condition. It is an upper bound of the same type as the 
lower bound in theorem 2.1, except that the quite natural rate 1/4 is replaced 
by the technical rate (resulting from the complex multiplier lemma 4.4): 

Theorem 2.2. For any a > a* defined by (4), there exists C > such that, 
for B = 1 or B = d s , for all L > 0, T G ]0, inf(vr, Lf] and u G L 2 (0,L), 
there is a g G L 2 (0,T) such that the solution u G C°([0, oo), L 2 (0, L)) of the 
following heat equation on [0, L] controlled by g from one end: 

d t u~d 2 s u = in ]0,T[x]0,L[, (Bu) ]s=Q = , U]s=L = g , u ]t =o = u , 

satisfies u = at t = T and ||g||L 2 (o,T) < Ce aL ^ t ||uo||l 2 (o,l) • 

Theorem 3.1 in [Sei84] yields this theorem for a* = 4/3* with (3 # 42.86. 
This result of Seidman can be improved to a* = 8(3* with /3* ~ 4.17 using 
his theorem 1 in [Sei86]. The value a* defined by (4) in theorem 2.2 is the 
best we obtained yet following the well trodden path of the harmonic analysis 
of this problem (cf. [Rus78] and [SAI00] for seminal and recent references). 
As explained at the end of the previous subsection, we conjecture that a* = 
1/4 is the optimal rate. The proof of theorem 2.1 also applies here, so that 
theorem 2.2 does not hold with a* < 1/4. This theorem is valid for more 
general linear parabolic equations and boundary conditions as formulated in 
theorem 4.1. 

2.3 Upper bound under the geodesies condition 

Our third result gives a good reason to strive for the best rate a* in theo- 
rem 2.2. In section 5, we prove that the upper bound for the null-controllability 
cost of the heat equation on a segment controlled at one end — the particular 
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case in which the computation are the most explicit — is also an upper bound 
for the multidimensional case of equation (1) under the following geodesies 
condition on the control region: every generalized geodesic in M intersects fl. 

In this context, the generalized geodesies are continuous trajectories t i— > x(t) 
in M which follow geodesic curves at unit speed in M (so that on these 
intervals 1 1— > x(t) is continuous); if they hit dM transversely at time to, then 
they reflect as light rays or billiard balls (and t i— > x(t) is discontinuous at 
to); if they hit dM tangentially then either there exists a geodesic in M which 
continues t h- > (x(t),x(t)) continuously and they branch onto it, or there is no 
such geodesic curve in M and then they glide at unit speed along the geodesic 
of dM which continues t h- > (x(t), x(t)) continuously until they may branch 
onto a geodesic in M. For this result and whenever generalized geodesies are 
mentionned, we make the additional assumptions that they can be uniquely 
continued at the boundary dM (as in [BLR92] , to ensure this, we may assume 
either that dM has no contacts of infinite order with its tangents, or that g 
and dM are real analytic), and that fl is open. 

Theorem 2.3. Let be the length of the longest generalized geodesic in M 
which does not intersect fl. If theorem 2.2 holds for some rate a* then the 
null- controllability cost of the heat equation for small time (cf. definition 1.1) 
satisfies the following geometric upper bound: 

lim sup T In Cr,n < a*L^ (5) 

When comparing this result to the lower bound in theorem 2.1, one should 
bear in mind that Ln is always greater than 2 sup^g^- dist(y, fl) (because the 
length of a generalized geodesic through y which does not intersect fl is always 
greater than 2 dist(y, fl)) and can be infinitely so. For instance, on the sphere 
M = S n , if fl is the complementary set of a tube of radius e around the 
equator, then sup ygM dist(?/, fl) — e and = oo. If fl is increased by a tube 
slice of small thickness 5, then the first length is unchanged while the second 
length becomes greater than the length of the equator of M minus 5, so that 
Lq is finite yet much greater than sup ygM dist(?/, fl) as e — > 0. 

Moreover, as recalled in section 1, this geodesies condition is by no means 
necessary for the null-controllability of the heat equation. It is more relevant 
to the wave equation on M, for which it is a sharp sufficient condition for 
exact controllability in time T by interior controls on fl as proved in [BLR92] 
(cf. theorem 5.3 for the precise statement). It was later proved in [BG97] that 
this condition is also necessary when the characteristic function of ]0, T[xfl is 
replaced by a smooth function 9 such that {9(t,x) ^ 0} =]0, T[xfl. 

In fact we use the exact controllability of the wave equation to prove our result 
on the null-controllability of the heat equation. This strategy was already 
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applied by Russell in 1973, but he used a complex analysis detour (cf. [Rus78]). 
In [Sei84], Seidman applied Russell's method to obtain an upper bound which, 
taking [BLR92] into account, corresponds to theorem 2.3 with a* = /5* ~ 
42.86. Theorem 2.3 improves Seidman's result beyond this slight improvement 
of the rate a* insofar as the complex analysis multiplier method he uses does 
not necessarily allow to reach the optimal a* in theorem 2.2. 

The control transmutation method (cf. [Her 75] for a survey on transmutations 
in other contexts) introduced in section 5 relates the null-controllability of 
the heat equation to the exact controllability of the wave equation in a direct 
way (as opposed to Russell's indirect complex analysis method). It is well- 
known that the geometry of small time asymptotics for the homogeneous heat 
semigroup t i— > e tA on L 2 (M) can be understood from the even homogeneous 
wave group t h- > W(t) (i.e. the group defined by W(t)wo = w(t) where w solves 
equation (53) with / = and Cauchy data (w, d t w) = (wq, 0) at t — 0) through 
Kannai's formula (cf. [Kan77], [CGT82], and section 6.2 in the book [Tay96]) : 



Our main idea is to replace the fundamental solution of the heat equation on 
the line e~ s2 l^ j^Tt appearing in Kannai's formula by some fundamental 
controlled solution of the heat equation on the segment [-L, L] controlled at 
both ends. We use the one dimensional theorem 2.2 to construct this funda- 
mental controlled solution in section 5. 

2.4 Open problems 

We shall now survey some questions raised by the results we have presented 
which we have been unable to answer yet. 

To improve the rate a* in theorem 2.2 by a complex analysis method, one 
could use the first method in [FR71], i.e. compute the null-controllability cost 
on the half-line [0, +00) explicitly by Vandermonde determinants and prove 
a quantitative version of Schwartz's theorem in [Sch43], i.e. estimate with 
respect to L the best constant cl in the following statement : every u in the 
closed linear hull in L 2 (0, +00) of the real exponential sums t 1— > e~ k t (k e N*) 
satisfies ||w|| L 2 (0i+oo) < c L ||u|| L 2 (0)i) . 

Theorem 2.3 opens new tracks to improve the upper bound for the null- 
controllability cost of (1) under the geodesies condition by methods which 
are not complex analytical. To improve the rate a* in theorem 2.2 (or in the 
multidimensional case of equation (1) when Q and M are star-shaped with re- 
spect to the same point) one could adapt the variational techniques (e.g. the 
log convexity method) or the Carleman's inequalities devised to prove unique 




.2 



/i4t) W{s) ds . 



(6) 
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continuation theorems. 



In the general case (without the geodesies condition), one could try to adapt 
the null-controllability proofs which use Carleman inequalities with phases 
to obtain an upper bound similar to the lower bound in theorem 2.1 in terms 
of the following distance function d : d(x, y) = sup{0(y) — <f>(x)}, for all x and 
y in M, where the supremum is taken over all Lipschitz functions : M — > K. 
with |V0| < 1 almost everywhere. There is a more geometric characterization 
of d in terms of path of least action (cf. section 2 of [Nor97]). 



3 Lower bound 

The purpose of this section is to prove theorem 2.1. 

As in section 1, let O be an open set in the n-dimensional Riemannian manifold 
M such that Q C M. Let (oUj)j e ^* be a nondecreasing sequence of nonnegative 
real numbers and (ej)j e ^* be an ortho normal basis of L 2 (M) such that e 3 - is 
an eigenvector of —A with eigenvalue u?. The heat kernel k can be defined for 

all t > and (x,y) G M by k(t,x,y) = J2j ex P(~toOj)ej(y)ej(x). Our main 
ingredient is Varadhan's formula which says that (cf. th. 1.1 in [Nor97] for 
example): 

\imt\nk(t,x,y) = —d(x,y) 2 /4 uniformly on compact sets of M 2 . (7) 

We shall also use Weyl's asymptotics for eigenvalues: 

3W > 0, #{j G N* | u)j <u}< Wu n (8) 

and the following consequence of Sobolev's embedding theorem: 

BE > 0, Vj G N*, He^lioo < Euf 2 (9) 

(cf. section 17.5 in [H6r85] for example). The unique continuation property for 
elliptic operators implies that Y = {y G M \ Q \ ei(y) ^ 0} is an open dense 
set in M \ Q, so that the supremun in theorem 2.1 can be taken over y G Y 
instead of y G M. 

Let y EY and a < d(y, Q) 2 /4 be fixed from now on. To prove theorem 2.1 we 
shall find A > and, for all T g]0, 1] small enough, some data G L 2 (M) 
such that ||e* A Uo||L 2 ((o,T)xn) — Ae~ a / T \\e TA UQ ||l 2 (jw)- To give further insight 
into the problem, we shall construct each Uq as a linear combination of a finite 
number of modes tj only. 

Let (3 be a real number such that a < j3 < d(y, f2) 2 /4. Since fix {y} is compact 
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in M 2 , Varadhan's formula (7) yields real numbers B > and T g]0, 1] such 
that 

V* G]0,T], Mx G H, \k(t,x,y)\ < Be'^ 1 . (10) 
Let e g]0, 1] small enough as specified later. For all T e]0,T/(l consider 
the data Uo(x) = Sw^eT)- 1 ex P( — £ Toj?j)ej(y)ej(x). To estimate the corre- 
sponding solution 



T 



(t, x) = (e tA <) (x) = X! exp(-(eT + t^e^e^x) , 

Wj<(eT)-i 



we compare it with k(eT + t,x,y). Using that the heat semigroup is a con- 
traction on L 2 (M), Parseval's identity and (9), we obtain 



sup \\k{eT + t,x,y) - u T (t,x)\\ L 2 (M) < \\k(eT, x, y) - Uq (x)\\ L 2 {m) 
te]o,T] 

= E \e~ £T ^ ej (y)\ 2 <E E e~^<E' ]T e^'/ 2 , 
for some £" > 0. But, Weyl's law (8) yields, for c > c > and 7 > 7o > 0, 



J2 e-T"* = E E e ~ 7 " J < W E ((^ + l)c) n e~ fec7 

Wj>c feGN* kc<ujj<(k+l)c feeN* 

< Vy 7Q ^ e -kc~f e (k+l)c~f/4 = W ^ e ~cj/2 J- e ~3kcy/i < 

where W l0 and W 7 ^^ are positive real numbers which depend on their indices 
but not on c and 7. Hence, with c = (eT)" 1 > 1 = Co and 7 = 70 = 1/2, we 
obtain: 

35' >0, VtE]0,T] \\k(eT + t,x,y) - u T (t,x)\\ L 2 {M) < B'e' 1 '^ 

Together with the estimate on k(eT + t,x,y) which follows from (10), this 
estimate yields by the triangle inequality, choosing e < 1/(4/3) and setting 
B" = \n\V 2 B + B', 



u 



T| |^((0,T)xQ) < (T|fi|) 1 /2 jBe -/3/((l+^) +T l/2 B / e -l/(4 £ T) < s „ e -^/((l +e )T) 



But using Parseval's identity and y G F, we have 

\\e TA u T J LHM) = [ £ le-^^e,^)! 2 ) > e^\ ei (y)\ > . 

Hence, choosing 5 small enough so that ct < (3/(1 + e) and setting 
A = e _2a; i|ei(y)|S", we have 



VT e]0,T/(l + e)], |K T |U 2(( o,T)xQ) < Ae- a / T ||e T Xll^ 



m ■ 
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Since A does not depend on T, this ends the proof of theorem 2.1. 



4 The segment controlled at one end 

In this section we prove theorem 2.2 for a more general linear parabolic equa- 
tion on a segment controlled at one end (in particular, it proves that the- 
orem 2.2 is true for the heat equation on a segment with any Riemannian 
metric). We follow [FR71] quite closely. 

For a positive a control time T, we consider the following mixed Dirichlet- 
Cauchy problem on the space segment [0, X]: 

d t u = d x (p(x)d x u) + q(x)u for (t, x) g]0, T[x]0, X[ , (11) 
(ao + b d x ) M lx=0 = , (ai + bid x ) n lx=x = g , u^ t=0 = u , (12) 
al + b 2 = a\ + b\ = 1, <pe C 2 {[0,X}), q G C°([0, X]) . (13) 

With assumptions (13), the operator A on L 2 (0, X) with domain D(A) defined 
by 

(Au)(x)=d x (p(x)d x u(x)) + q(x)u(x) 

D(A) = H 2 (0, X) n {(a + b d x ) u ]x=0 = (a t + b±d x ) u ]x=x = 0} 

is self-adjoint and has a sequence {— A n }„ e N* of increasing eigenvalues and an 
orthonormal Hilbert basis {e n } ne ^* in L 2 (0,X) of corresponding eigenfunc- 
tions, i.e. : 

Vn G N*, -Ae n = \ n e n and A n < X n+1 . 
Moreover, (13) ensures the following eigenvalues asymptotics (cf. [FR71]): 

3u G K, X n = — (n + v) + 0(1) as n — > oo , where L — / yp(x) dx . 
L Jo 

(14) 

Theorem 4.1. For any a > a* defined by (4), £nere exists C > sncn £/ia£, 
/or any coefficients (13), /or a// T G ]0, inf (ir, L) 2 ] and uo G L 2 (0,X) there is 
a control g G L 2 (0,T) such that the solution u G C°([0, oo), L 2 (0, X)) o/ (11) 
and (12) satisfies u = at t = T and || f? || z, 2 (o,^r) < Ce aL ^ T ||no||L 2 (o,x) • 

As in [FR71], the proof applies to the slightly more general eigenvalue asymp- 
totics A n = ^r- (n + i/) + o(n). We divide the proof of this theorem in three 
steps. 
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4-1 Reduction to positive eigenvalues, to a segment of p-length L = ix , and 
to the control window ] — T/2, T/2[ 

As a first step, we reduce the problem to the case Ai > by the multiplier 
t i— > exp(At), to the case L = n by the time rescaling 1 i— > at with a = (tt/L) 2 , 
and to the time interval [— T/2, T/2] by the time translation t i— > t — T/2. 

The function u satisfies d t u = Au and (ai + u-\ x =x = g if and only if 
u(t,x) = exp(Xt)u(t,x) satisfies d t u = Au and (a± + bid x ) u-\ x =x = 9 with 
A = A + A and = exp(Xt)g(t). For any A > — Ai, the lowest eigenvalue 
of A > Ai + A > is positive. In A, q is changed into q + A and p is un- 
changed so that L is unchanged. Moreover ||<7||l 2 (o,t) < ex P(^/2)||^||L 2 (o,T) so 
that ||^||l 2 (o,t) ^ C' eaL2/ " r ||' u o||L 2 (o,x) implies the estimate in theorem 4.1 with 
C = Cexp(A7r/2). This proves the reduction to positive eigenvalues. 

We now prove the second reduction. Assume the theorem is true when L takes 
the value L = ir. Given L > and T G ]0, inf (ir, L) 2 ] we set f = a 2 T e]0, L 2 ] 
and A = o" 2 A, where a = (tt/L) 2 . By applying the theorem to A on ]0,T[, we 
obtain H^H^^t) — Ce aL2 ^ T \\uo\\L2^ 0tX )- The function g(t) = g(at) is a control 
for the solution u(t, x) = u(at, x) of d t u = Au on ]0, T[ at the cost ||g||L 2 (o,T) = 
||<?IU 2 (o,T)-k/ 7r - Since T < ir 2 implies L/n < (L 2 /T) 1 / 2 , for all a > a there is 
a C such that for all L > and T G ]0, inf(7r, L) 2 ]: Ce^'^L/ir < Ce aL ^ T . 
Therefore g satisfies the estimate in theorem 4.1. 

These two reductions allow us to assume from now on Ai > and L = n. 
Making a weaker assumption on the remainder term in (14), we shall only use 
the following spectral assumption: 

Vn G N*, < A„ < A n+1 and 3v G R, A n = (n + v) 2 + o(n) as n — > oo . 

(15) 

It is obvious that theorem 4.1 is invariant by time translations and we shall 
prove it for the control window ] — T/2, T/2 [ instead of ]0, T[. 

4-2 Spectral reduction to a problem in complex analysis 

In this second step, we recall that the control g in this theorem can be ob- 
tained as a series expansion into a Riesz sequence {g n }nm* m L 2 (—T/2,T/2) 
which is bi-orthogonal to the sequence {exp(— A n t)} neN ». We also recall how 
the Paley- Wiener theorem reduces the construction of such biorthogonal func- 
tions to the construction of entire functions with zeros and growth conditions 
(this well-known method in complex analysis is the second method in [FR71] 
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called the Fourier transform method there). Our estimate on the control cost 
1 1 <7 1 1 l 2 f-T/2,T/2) relies on a good estimate of \\g n \\L 2 (-T/2,T/2) as T tends to zero. 
This additional difficulty was first taken care of by Seidman in [Sei86] for 
A n = in 2 and it was recently overcome for more general sequences in [SAI00]. 
Our contribution is a slight improvement on the estimates of Seidman and his 
collaborators in our less general setting. 

In terms of the coordinates c = (cj)j e ^* of u$ in the Hilbert basis (ej)j £ ^*, the 
controllability problem in theorem 4.1 is equivalent to the following moment 
problem (by straightforward integration by parts, cf. [FR71]): 

e-^-SnSW dt = -e-^ T c n , 

J-T/2 

where 7ri = e n (X)p(X)/b 1 if 61 ^ and 7ri = -e' n (X)p(X)/a 1 if b± = 0. In 
both cases, the asymptotic expansion of e n yields that (| 7n |) is bounded from 
below by some positive constant 7. If {g n }neN* i n L 2 (—T/2, T/2) is a sequence 
which is bi-orthogonal to the sequence {exp(— A n t)} n£ N*, i.e. 



9n{t)e- Xnt dt = 1 and VA; G N*, k ^ n, / ^„(t)e" Afc * dt = 0, 

-T/2 J-T/2 

(16) 



00 c 

then g(t) = — —e~ XnT ^ 2 g n (—t) is a formal solution to this moment prob- 

n=l 7n 

lem. The following theorem in complex analysis allows to construct a bi- 
orthogonal sequence such that this series converges and yields a good estimate 
of \\g\\L 2 (~T/2,T/2) as T tends to zero. 

Theorem 4.2. Let a* be defined by (4). Let {A n } n£ N* be a sequence of real 
numbers satisfying (15). For all e > there is a C £ > such that, for all 
t e]0, 1] and nGf, there is an entire function G n satisfying 



G n is of exponential type r, i.e. limsupr 1 sup In |Gr n (;z)| < r, (17) 

\z\=r 

G n {i\ n ) = 1 and \/k eW, n, G n (i\ k ) = 0, (18) 

(p-\-oo \ 1/2 

j \G n (x)\ 2 dx) <C £ e £ ^e a *^ +2£ > 2 /^ (19) 



According to the Paley- Wiener theorem (1934), (17) implies that the function 
x 1— > G n (x) is the unitary Fourier transform of a function t 1— > g n (t) in L 2 (IR) 
supported in [— r, r]. With r = T/2, this yields: 

1 r T / 2 

G n (x) = ^= g n {t)e- ltx dt and ||<y i2 = ||G n || L2 . (20) 

V27T J-T/2 
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Hence (18) implies (16) and (19) implies that the series defining g converges 
with: 



l#IU 2 < 



1/2 



e- X - T ' 2 \\g n \y < \\u4vi-e a ^ +£)2/T £ 



n=l 

Since as T — > we have 
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Vn=l 



£ e - A " r e 2£ v^ < e ^/T e ~ x - T / 2 ~ e 2e2 / T (T/2)- 1 / 2 r(l/2)/2 < C> 3£2 / T , 



n=l 



n=l 



this implies ||fi'||L 2 (-T/2,T/2) < Cae Q7r2/T ||wo|U 2 (o,x), with a = + 2e/7r) 2 + 
3£ 2 /7r 2 and C a = C £ C' £ /^. Since ct — > a* as £ — > 0, this completes the proof 
that theorem 4.2 implies theorem 4.1. 



^.5 Complex analysis multipliers 



In this subsection we shall prove theorem 4.2 by the following classical method 
in complex analysis (cf. section 14 in [Red77] for a concise account with ref- 
erences, and the two volumes [Koo92] for an extensive monograph on mul- 
tipliers): for all n G N* and small r > 0, we shall form an infinite product 
F n normalized by F n {i\ n ) = 1 with zeros at i\k for every positive integer 
k 7^ n, and construct a multiplier M n of exponential type r with fast decay 
at infinity on the real axis so that G n = M n F n is in L 2 on the real axis. At 
infinity, it is well known that the growth of z i— > ln|F n (z)| can be bounded 
from above by a power of \z\ which is inverse to that of n \— > \i\ n \ ~ n 2 (cf. 
theorem 2.9.5 in [Boa54]) we prove that our lnF n is essentially bounded by 

z I— > 7TyJ\z\ + o(y/X^) where the constant n is optimal (cf. remark 4.5). There- 
fore M n has to be essentially bounded by C n {r) exp(— it\J\x\) on the real axis, 
for some constant C n (r) > 0. The key point (as in [Sei84], theorem 1 in [Sei86] 
and theorem 2 in [SAI00]) is to construct a multiplier M n such that C n (r) has 
the smallest growth as r tends to 0. The following two lemmas give the key 
to the construction of F n and M n respectively. 

Lemma 4.3. Let {A n }„ GN » be a sequence of real numbers satisfying (15). For 
all e > there is a A £ > such that, for all n G N* ; the entire function f n 

defined by f n (z) — ] [ ( 1 — — ) satisfies 



ln|/n(*)l < (n + e)^\z\+A e (21) 
|ln|/n(A n )|| <eJ\ n + A e (22) 
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Proof. For every n £ N*, we introduce the counting function of the sequence 

{^k}k<=N*\{n}'- 

N n {r) = #{k G N* \ {n} | \ k < r} . 

From (15) we have N — 1 < N n < N and y/X^ = n + z/ + o(l). Since 
Afe < r < A fc+ i implies - k < yfr - N (r) < \/X k +i - (k + 1) + 1, we 
deduce | y/r — N n (r) — u\ < 2+o(l). The proof uses the assumption (15) through 
the estimates of the increments A n := \ n+ i — X n and A n := y/\ n +i — \fKi and 
their increments: 



A n = n 2 + 2un + o(n) , A n = 2n + o(n) , A n - A n _i = o(n), 

= n + v + o(l) , A n = l + o(l), A n - A n _ x = o(l) , 
VrG]0,A 1 [,iV n (r) = 0, > 0, Vr, jy^ - N n (r)\ < A . 

We shall use repeatedly that for any real sequence {r n } TteN . 



r n = o(l) 



In 1 + 



k»|(i + o(i)) 



1 + (1), 

To prove (21), we estimate the left hand side in terms of N n : 



(23) 
(24) 
(25) 



(26) 



ln|/ B (*)|<5> l + Y 1 

k^n \ Ak / 



= / Njr 
Jo 





z 


dr 


\z 


+ r r 



In \l + L±\ d N n {r) 



N n (\z\s) ds 



1 + s s 

To estimate this last integral we use (25) and the integral computations: 



r°° y/s ds _ r°° 2dr r 
Jo l + s~ ~ Jo 1 + r 2 ~ 77 ' 7? 



00 2dr 



r °° ds 

I ft 8(1 + 8) 



In 


s 




l + s\ 



*1 

FT 



ln(l + 



Thus we obtain In \ f n (z)\ < nJ\z\ + Aln(l + J^), so that, for all e > there 



is a A' £ > such that In \f n (z)\ < (n + e)^/ |z| + 

To prove (22), we estimate the left hand side in terms of N n : 



ln|/„(A„)| = £ln(^-l)+£ln ( 1 - ^ 

A k J \ A fe , 



k<n \ fc / k>n \ 

K -' In - l) dN n (r) + r In { 1 - ^) diV„(r) 



n+l 
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Integrating by parts yields In |/ n (A n )| = I n + B n with 



4=/ Wn(r)T— = +/ iV n r 

A r Xn-r r Jx- +1 



X„ dr f 00 



A ra dr 
X n — r r 



Br, 



N n {r)\n 



+ 



7V n (r) ln(l-y 



l n+l 



To estimate the boundary term S n , we first simplify its expression using 
N n (Xi) = and N n (Xn-i) = N n (X~ +l ) = n — 1, then we sort out the in- 
crements A n = A n+ i — X n , and finally we use (23) and (26): 



B n = (n-1) 
= (n-l) 
= (n-l) 



In 



A, 



, A n _i 



1 -In 1 



A r 



A 



n+l , 



ln|l-^Uln(l + A " + Vl ' 



A, 



A n -i 



^(i + o(D) + *i±42>(i + o(D) 



= o(D . 



Now we estimate the integral term /„. Performing the change of variable r 
X n s and using (25) yields: \I n — \/X^J n \ < AK n with 



(is 



Jn 4" (l-^)v^ + ^(l- S )v^ 



(is 



(l-s)s y-S±i(s-l)s 



(is 



The term is readily computed and estimated using (23): 



In 



2 In 



1 - s 

An 

n 2 + 0(n) 
2n + o(n) 



+ 



In 



s - 1 



^= ln ¥ +ln ^ +ln ( A -4 + 

— H — - il n il -n— 1 \ /x l 



+ 21 nv /A n + 0(l) = o(JA n ) 



We compute J n after a change of variable, and estimate it by (24) and (26) 
after sorting out the increments A n = a/A„ + i — V%i : 
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J, 



2dr 
r 2 — 1 



2dr 



In 



1 — r 



In 



+ VX 



-In 



r + 1 



+ 



In 



r — 1 
r + 1 



V A n + l 



(l- A » A A -'Ulnfl + - A " + A 



In 1 



2VAi 



n-l , 



A 

o(l)(l + 0(1)) + ^11(1 + o(l)) + ^(1 + 0(1)) = 0(1) . 



2n 



n 



Plugging the estimates K n = o(y/X^) and J n = o(l) into \I n — \f\^J n \ < AK n 
yields I n = o{\f\ l ). Plugging this estimate and B n = o(l) into In |/ n (A n )| = 
I n + B n yields In |/„(A„)| = o( v / A^), which completes the proof of (22). □ 

Lemma 4.4. Let a* be defined by (4). For all d > there is a D > such 
that for all r > 0, there is an even entire function M of exponential type r 
satisfying: M(0) = 1 and 

a d 2 

Wx>0, ln\M(x)\ < ^— + D - dVx and \M(ix)\ > 1 . (27) 

2r 



Proof. Following Ingham and many others since 1934 (cf. section 14 in [Red77] 
for theorems and references) we seek a multiplier M of small exponential type 
decaying rapidly along the real axis in the following form: 



M(z) = Y[ smc ( ~~~) where sinc(O) = 1, Vz e C*, sinc(,z) 



sin(z) 



(28) 



and where {a„} nG N is a non decreasing sequence of positive real numbers such 
that tm = SneN ~t~ < 00 • Since the cardinal sine function sine is an even 
entire function of exponential type 1 satisfying sinc(O) = 1 and sinc(ix) = 
sinh(x)/a; > 1 for all x > 0, (28) defines an even entire function M of expo- 
nential type r M satisfying M(0) = 1 and |M(ia;)| > 1 for all x > 0. 

We define {a n }„ G N by the slope A of its counting function N and its first term 
a (to be chosen large enough): 

N(r) := Yl 1 = i A Vu\ for r > 2 and a > A~ 2 , 

\a n \<r 

where [x] denotes as usual the greatest integer smaller or equal to the real 
number x. The exponential type tm of M is easily bounded from above by 
r = 2A/^: 



™ := £ f 



diV(r) 



7V(r) 



o r z 



, 2A 
— — <ir = =: t 
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and we are left with estimating the decay of: 

In \M(x) | = r f (-) dN{r) where f(9) = In sinc(#) = In . (29) 

We shall choose A such that, for all a > A~ 2 , In |M(x)| < —d^fx + 0(1) as 
x — > +00, and then prove that: In |M(x)| < a*d 2 /(2r) — rf v /a : + 0(l) as r — > 
(equivalently a — > +00) uniformly in x > 0. 

For x > a we take advantage of the boundedness of sine through the estimate 
f{9) < — In \9\ for \9\ < 1, by splitting the integral in (29) into the two terms: 



'=£-'(7) dN(r)i M dN(r) -£ N(r) T 

J = J™f(f) dN ( r ) = I'/'W^ (I) <0-/(iW*) 

where right hand sides were integrated by parts and 9 = x/r. Now we plug in 
the basic estimate on AT: A^fr — 1 < iV(r) < A^fr for r > 2. The first term is 
now estimated by 

rx rx fly , . 

I<-A -7=+/ — = -2A(y/x- y/aoj +lna;-lna . (30) 

Jao yT Jao T ^ ' 

To estimate the second term, we first observe that the Hadamard factoriza- 
tion of the cardinal sine function sinc(7r,2) = IlneN* (l ~~ ^s) an d the Taylor 
expansion of the logarithm at 1 imply: 

m = _j2 iffi (T^ for |0| < 1 , where = E ^ ■ 

The second term is now estimated by 



j<£f'{e)(^£-i\de-f(i)AjZ (31) 
= -^E (5^-1)^-/(1) = -^-^- /(i), 

where the series for / was differentiated, multiplied and integrated term by 
term, and S* = EfeeN* k (4k-i) C ffi - Putting (30) and (31) together yields: 

Vx>a , ln|M(x)| < -(2 + T?)Ay/x + lux - /(l) + 2A^ . 
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so that, for all d > (2 + T,*)A there is a D 1 such that 
\/d > (2 + S*)A, 3D 1 >0,\/x> oq, In |M(x) | < 2A^ - + D 1 . (32) 



Since | sine | is bounded by 1: for all x, In |M(x)| < 0. Moreover d > 2 A, so 
that (32) implies 

Va > A~ 2 ,Vx > 0, \n\M(x)\ < d^/E~ Q - d^fx~ + D x . (33) 

Since d > (2 + S*)A and r = 2A/ v / a ~, this proves: 

Vr < 2A 2 ,Vx > 0, \n\M(x)\<^--d^ + D 1 (34) 

2r 

with ai = 4/(2 + 11*). 

For x < a , we can also use the better estimate: 

In |M(x)| < y / (-J rfiV(r) = jf /'(0)7V (-J (35) 

k(4k - 1) \a irj 

If we keep only the first term (i.e. k — 1) of the series in (32) and (35), we get 
that for all d > (2 + \^r)A there is a L> 2 such that : 

\/x> a , ln\M{x)\<2A v ^-d^ + D 2 (36) 
\/x < a,,, In |M(x)| < -A^^ - / (1) . 

Now, for all x < a 



In |M(x) | - 2A^d~Q + < Ay/d^F (—) 

Veto/ 



withF(X) = -2+(2+e)VX+l^(VX-AX 2 ) = -2+(37/18+e) v / X-2X 2 /9 
and e = d/A— (2+|^) > 0. Since F is increasing on [0, 1] and F(l) — e—1/6, 
choosing A so that e < 1/6, yields that In |M(x)| — 2A v /a " + dy^ < 0, for all 
x < a . Together with (36), this proves 
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Vx > 0, In \M{x) | < 2A^/E~ - d^ + D 2 (37) 
Since d > (2 + \^)A = 37A/18 and r = 2A/ V ^, this proves: 

Vr < 2A 2 ,Vx > 0, In | M(x) | < - dy/x + D 2 (38) 

2r 

with a 2 = 2(36/37) 2 . 

Equations (34) and (38) complete the proof of the lemma 4.4 with a* = 
minjai, a 2 }. Since we have checked on a computer that ot\ > a 2 , we decided 
to state the lemma with a* = a 2 , i.e. (4). □ 

To prove theorem 4.2, we use lemmas 4.3 and 4.4 with d = it + 2e and define : 
G n = F n M n with F n (z)=f n (-iz)/f n (\ n ) and M n (z) = M(z)/M(i\ n ) . 

Thanks to lemma 4.3, the entire function F n satisfies 

F n (i\ n ) = 1 and V/c G N*, k ^ n, F n (i\ k ) = 0, (39) 
\n\F n (z)\<(n + e)^/\z\+eJ\ l + 2A e (40) 

where (39) is an obvious consequence of the definitions of f n and F n , and (40) 
is a consequence of the estimates (21) and (22). 

Thanks to lemma 4.4, there is a D £ > such that the entire function M n is 
of exponential type r and satisfies 



M n (i\ n ) = 1 (41) 

ad 2 i — 
ViGR, \n\M n (x)\ < -^— + D £ -dJ\x\ (42) 

2t v 

where (41) is an obvious consequence of the definitions of M and M n , and 
(42) is a consequence of (27) since M is even. 

The entire function G n has the same exponential type as M n since (40) implies 
that the exponential type of F n is 0. Hence (17) holds. Putting (39) and (41) 
together yields (18). Since d = it + 2e, (40) and (42) imply 

/ — / — ad 2 
VieR, \n\GJx)\ < D £ + 2A £ -eJ\x\ +eJ\ n + 



is 



2r 

»+oo f— \ !/ 2 



Hence (19) holds with C £ = e D " +2A " ( / e" 2e v 1*1 dx J . Theorem 4.2 
proved. 
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Remark 4.5. Under assumption (15), lemma 3 in [SAI00] (which applies to 



(40), the estimate 0{\/\,) improves to o(y/X^) and the constant 2tt improves 
to the optimal ir (optimality can be deduced from theorem 4.1.1 in [Boa54]). 

Seidman obtained lemma 4.4 for a* = /3* with /3* ~ 42.86 in the proof of 
Theorem 3.1 in [Sei84]. His later Theorem 1 in [Sei86] improves the rate to 
a* = 2/3* with /3* ~ 4.17. Theorem 2 in [SAI00], which applies to much more 
general spectral sequences, yields lemma 4.4 for a* = 24. The argument used 
in section 3 can be used to prove that lemma 4.4 does not hold for a* < 1/4. It 
would be interesting to determine the smallest value of a* for which it holds. 



5 Upper bound under the geodesies condition 

In this section we prove theorem 2.3 in three steps. T>'(0) denotes the space of 
distributions on the open set O endowed with the weak topology and A4(0) 
denotes the subspace of Radon measures on O. 

5.1 The segment controlled at both ends 

In a first step we prove that the upper bound for the null-controllability cost 
of the heat equation on the segment [0, L] controlled at one end is the same 
as the null-controllability cost of the heat equation on the twofold segment 
[-L, L] controlled at both ends. 

Given a time T > and a length L > 0, we denote by D (respectively N) 
some continuous operator from L 2 (0,L) to L 2 (0,T) allowing to control to 
zero in time T the heat equation on [0, L\ with zero Dirichlet (respectively 
Neumann) condition at by a Dirichlet control at L. More precisely, for all 
u G L 2 (0,L) the solution u G C°([0, oo), L 2 (0, L)), denoted by u = SdUq 
(respectively u = SnUq), of the Cauchy problem in theorem 2.2 with B = 1 
(respectively B = d s ) and g = Duq (respectively g = Nu ) satisfies u — at 



Proposition 5.1. For any time T > and any length L > 0, there is a 
continuous operator K from L 2 (—L,L) to L 2 (0,T) 2 allowing to control to 
zero in time T the heat equation on [— L, L] by Dirichlet controls at both 
ends at the same cost as D and N, i.e. for all t> G L 2 (—L,L) the solution 




t = T. 
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veC°([0,oo),L 2 (-L,L)) of: 
d t v-d 2 v = in]0,T[x]-L,L[, (v ]s= . L , v ]s=L ) = Kv , v ]t=0 = v (43) 
satisfies v = at t = T and \\K\\ < s\xp(\\D\\, \\N\\). 

Proof. Given v G L 2 (—L,L), we decompose it in odd and even parts : 
= v 0tOdd + v 0teven . We denote by u QMd and u^ even the restrictions of v 0tOdd 
and fo,ei)en to [0, L], We denote by / = Duo t0dd and g = NuQ t£ven the corre- 
sponding controls. We denote by u odd = S D u 0>odd and u even = S N u , e ven the 
corresponding solutions. 

We define v G L 2 ([0,T] x [-L, L]) by v(t,±s) = u even (t,s) ±u odd (t,s) for 
s > 0. Since 

(d t - dl)u even = (d t - d 2 )u odd = in V'(]0,T[x]0,L[) , 

we have, denoting the Dirac mass at s = by 5 S G X>'(R), 

- <9> = 2u odd (t, 0) <g> £(0) + 2d s u even (t, 0) <g> <f,(0) . 

But u odd (t,0) = d s u even (t,0) = by the definition of D and N. Hence 
(dt - d 2 s )v = 0. Moreover v(0,s) = v (s), v(T,s) = 0, v(t, L) = g(t) + f(t), 
v(t,—L) = g(t) — f{t). Therefore, setting Kv = (g — f,g + f) yields an 
operator K satisfying the null-controllability property required. 

To finish the proof we estimate its cost Taking the Euclidean norm for 

Kv G = (g-f,g + f),we have ||^o||!2 (0 ,r)2 = 2 \\f\\h(o,T) + 2 IMIl2(o,t)- Since 
/ = Du 0tOdd and g = Nu 0jeven , setting C = sup(||Z?||, ||iV||) we have 



ll^ollia^T) 2 <2C 2 (||uo,o<w||l2( 0)£ ) + ||wo,euen||l, 2 (0,L)J (44) 

Moreover, since 2u 0>odd (s) = v (s) - v (-s) and 2u 0)£ven (s) = v (s) + v (-s) 
for s G [0, L], we have 



||2«o,«id||i2 ( o,i) = \\vo\\ 2 L 2 hLjL) -2 v (s)v (-s) als (45) 

J 
r-L 

pu 0:even \\ 2 L 2 {m = \\v \\ 2 L 2 { _ LjL) + 2 J v (s)v (-s) als . (46) 
Equations (44), (45) and (46) imply \\Kv \\ L 2 ioX) 2 < C||v |U 2 (-i,i)- □ 
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5.2 The fundamental controlled solution 



In a second step we construct a "fundamental controlled solution" v of the 
heat equation on the segment controlled by Dirichlet conditions at both ends. 

Proposition 5.2. If theorem 2.2 holds for some rate a* ; then for any a > a*, 

there exists A > such that for all L > and T G ]0, inf(7r/2, L) 2 ] there is a 
v G C°([0, T],M(] — L, L[)) satisfying 

d t v - d 2 v = in £>'(]0, T[x] — L, L[) , (47) 
f]t=o = 5 and f]t =T = , (48) 

\MLH]0,T[x]-L,L[)< AeaL2/T - ( 49 ) 

We shall sometimes refer to a function v satisfying the above requirements as 
a fundamental controlled solution on ]0, T[x] — L,L\ at cost (A, a). 

Proof. We first reduce the problem to the case L = ir/2 using the rescal- 
ing (t,s) h-> (o- 2 t,as), cr > with a = n/(2L). Given L > and T G 
]0,inf(7r/2,L) 2 ], we set Z = 7r/2jmd f =^a 2 T e]0,L 2 ]. Let 6 be a funda- 
mental controlled solution on ]0, T[x] — L, L[ at cost (A, a). Setting v(t, s) = 
av(a 2 t l as) defines a fundamental controlled solution v on ]0, T[x] — L,L\ at 
cost (A/y/a,a). Since T < L 2 , we have A/y/a < A{L 2 /T) 1 / 4 . Hence for all 
a > a there is an A > such that v is also a fundamental controlled solution 
on ]0, T[x] — L,L[ at cost (A, a). Therefore, it is enough to prove proposi- 
tion 5.2 in the particular case L = ir/2. 

We assume theorem 2.2 holds for some rate a*. Let a > a* > a*, L = L = n/2 
and T g]0, L 2 ] be fixed from now on. We set a = (1 — e)<S* and T = (1 — e)T 
where e g]0, 1[ is chosen such that a > a*. Applying theorem 2.2 once with 
B = 1 and once with B = d s , and then applying proposition 5.1 yields a C > 
independent of T such that: 

\\K\\ < sup(||i}||, ||iV||) < Ce aL2/T = Ce & * L2/f . (50) 

We define v G C°([0, f], M{] - L, L[)) as the solution of 

d t v - d 2 s v = in ]0, f [x] - L, L[, (v ]s= _ L , v ]s=i ) = b, v ]t=0 = 5 

where the control b G L 2 (0,T) 2 is defined by b(t) = for t < ef and 
by b{eT + If) = K(v ]t=£T )(t') for t' G]0,T[. Note that v = v ]t=eT is just 
the Dirac mass at the origin smoothed out by the homogeneous heat semi- 
group during a time eT, so that Vq G L 2 (—L, L). Moreover sT + T = T and 
v(t, s) = v(eT + t, s) is the solution of (43), so that v^ t= f = v^ t =T = 0. 
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To finish the proof that v is a fundamental controlled solution on 
]0,T[x] - L,L[, we estimate its L 2 (]0,T[x] - L,L[) norm which we abbre- 
viate as \\v\\f i- Setting ej(s) = sin(j(s + 7r/2))y / 2/ 7r defines an orthonormal 
basis (ej)j e ^* of L 2 (] — L, L[) such that tj is an eigenvector of — A s with eigen- 
value j 2 . In the weak topology, the Dirac mass can be decomposed in this 
basis as S(s) = J2j e j(®) e j( s )- Note that the sequence (e J (0)) : , g N* is bounded. 
For t e]0, T], we introduce the coordinates (vj(t))j e ®* of v(t, ■) G L 2 (] — L, L[) 
in the Hilbert basis (ej)j eN ». Using these coordinates and abbreviating the 
L 2 (]0,T[) norm as || ■ \\f, the function v and its norm write 

rf 

V(t,s) = E fi j(*) e i( s ) alld NlfL=/ Sl^(*)| 2d * = Sll^llf ■ ( 51 ) 
3 ' J J 

As in [FR71], these coordinates can be computed by Vj(0) = e,(0) and 
^.(t) = e-^v^O) + jf e^*"* 7 ) ( e ;(-L)u(*', —L) - e'^LW , L)) dt' . (52) 

Using Young's inequality to estimate the second term of the right hand side, 
we have (since f < 4, |eJ(±L)| = \vj(0)\ = ^2/^ < 1) 

Nlf < l^(0)|||e-^|| f + \\e- j2t \\ Ll{]0>f[) (|e;.(-L)|||^,-L)|| f +|e;.(L)|||^',L)|| f ) 
<j(l + \\v(t',-L)\\ f +\\v(t',L)\\ f ) . 

Hence equation (51) implies 



11*% < (l + \W, -m + \\v(t', L)\\l) E J, = (l + ll^oll^Qcf D ) 

3 •> 

But there is an A' > independent of eT < 1 such that: 

ikoiii 2(] _ U[) =Ei^(^)r<E^' 2£f < A ' 



Hence equation (50) yields a C" > independent of T such that: 
WI W <^(1 + 2^11^11 hok)<^e^. 



Since a > a*, there is an A > independent of T such that: 
ll^llf i < Ae aL2,/T . This completes the proof that v is a fundamental con- 
trolled solution on ]0, T[x] - L, L[ at cost (A, a). □ 
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5. 3 The transmutation of waves into heat 



In a third step we perform a transmutation of an exact control for the wave 
equation into a null-control for the heat equation. Our transmutation formula 
can be regarded as the analogue of Kannai's formula (6) where the kernel 
e~ s Z( 4 *) /y/Airt, which is the fundamental solution of the heat equation on the 
line, is replaced by the fundamental controlled solution that we have con- 
structed in the previous step. To ensure existence of an exact control for the 
wave equation we use the geodesies condition of Bardos-Lebeau-Rauch (al- 
ready mentioned above theorem 2.3): 

Theorem 5.3 ([BLR92]). If L > L n then for all (w ,wi) G H^(M)xL 2 (M) 
and all (w 2 , w 3 ) G Hq(M) x L 2 (M) there is a control function f G L 2 (R + x M) 
such that the solution w G C°(R + , Hq(M)) n C 1 ^, L 2 (M)) of the mixed 
Dirichlet-Cauchy problem (n.b. the time variable is denoted by s here): 

d 2 w -Aw = l ]0 ,L[xn/ in R+ x M, w = on R+ x dM, (53) 

with Cauchy data (w,d s w) = (w ,Wi) at s = 0, satisfies (w,d s w) = (w 2 ,w 3 ) 
at s = L. Moreover, the operator S w : (H£(M) x L 2 (M)f -> L 2 (R + x M) 
defined by Sw ((wo,wi), (w2,ws)) = f is continuous. 

We assume that theorem 2.2 holds for some rate a*. Let a > a*, 
T g]0, inf(l, Lq)[ and L > be fixed from now on. Let A > and 
v G L 2 (]0,T[x] — L,L[) be the corresponding constant and fundamental con- 
trolled solution given by proposition 5.2. We define v G L 2 (IR 2 ) as the extension 
of v by zero, i.e. v(t, s) = v(t, s) on ]0, T[x] — L, L[ and v is zero everywhere 
else. It inherits from v the following properties 

d t v - d 2 v = in V'QO, +oo[x] - L, L[) , (54) 
v-\t=o = $ and V.-\ t =T = , (55) 
IMUwoo[xM) ^ ^ aL2/T • (56) 

Let Mo G Hq(M) be an initial data for the heat equation (1). Let w and 
/ be the corresponding solution and control function for the wave equation 
obtained by applying theorem 5.3 with wq = uq and ui\ = w 2 = w% = 0. We 
define w G L 2 (R; Hq(M)) and / G L 2 (R x M) as the extensions of w and 
/ by reflection with respect to s = 0, i.e. w(s,x) = w(s,x) = w(—s,x) and 
f(s,x) = f(s,x) = f(—s,x) on R + x M. Since w 1 = 0, equation (53) imply 

d 2 w - Aw = l]- LtL[x nf in V{R x M), w = on R x dM, (57) 
The main idea of our proof is to use v as a kernel to transmute w and / into 
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a solution u and a control g for (1). Since v G L 2 (R 2 ), w; G L 2 (R; Hq(M)) and 
/ G I 2 (l x M), the transmutation formulas 

«(t, X )= fv(t,s)w(s, X )ds and 9 (t,«) = / , (58) 

define functions u G L 2 (R;i?^(xM)) and # G L 2 (R x M). Since w(s,z) = 
d s w(s,x) = for \s\ = L, equations (57) and (54) imply 

d t u-Au= 1 g in£>'(]0,+oo[xM) and u = on ]0, TlxdM, (59) 
]o,T[xn 

The property (55) of v implies 

M]t=o = Mo and u-] t=T = . (60) 

Setting C = \^2A\\Sw\\, Cauchy-Schwarz inequality with respect to s, the 
estimate (56) and ||/||| 2(MxM) = 2\\S W ({u , 0), (0, 0))) ||| 2(R+xM) imply 

IblU^MxM) < ||^||l 2 (i; 2 )||/||l 2 (RxM) < Ce aL /T ||m ||j/1(m) • (61) 

We have proved that for all a > a* there is a C > such that for all 
m G Hq(M), T g]0, min{l, and L > Lq, there is a control c/ which 
solves the null-controllability problem (59), (60), at a cost so estimated in 
(61). The same property holds for the space of data L 2 (M) instead of Hq(M), 
since ||e £TA w ||^i (M) < \\u \\ L 2 {M) C / VsT with e g]0, 1[ and C = ||(1 + 
A)e _2v ^||^/ ]K - ) . Therefore limsupTlnCVQ < aL 2 . Letting a and L tend re- 
spectively to a* and Lq in this estimate completes the proof of (5). 
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